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Summary. — The quenching of jets (particles with pr » T,Aqcd) in ultra- 
relativistic heavy-ion collisions has been one of the main prediction and discovery at 
RHIC. We have studied, by a simple jet quenching modeling, the correlation between 
different observables like the nuclear modification factor Raa(pt), the elliptic flow 
V2 and the ratio of quark to gluon suppression R,AA(quark) / R,AA{gluon) . We show 
that the relation among these observables is strongly affected by the temperature 
dependence of the energy loss. In particular the large V2 and and the nearly equal 
Raa(j>t) of quarks and gluons can be accounted for only if the energy loss occurs 
mainly around the temperature T c and the flavour conversion is significant. Finally 
we point out that the efficency in the conversion of the space eccentricity into the 
momentum one (vz) results to be quite smaller respect to the one coming from 
elastic scatterings in a fluid with a viscosity to entropy density ratio Ann/ s — 1. 

PACS 12.38. Mh, 25.75 -q , 25.75 Ld - . 



1. Introduction 

The experiments at the Relativistic Heavy Ion Collider (RHIC) have given clear in- 
dications of the formation of Quark gluon plasma (QGP). One way to probe the QGP is 
to exploit the high energy jets (pr >> T,Aqcd) produced by the hard collisions at the 
initial stage. They are internal probes propagating through the fireball and interacting 
with the medium losing energy, hence carrying information on its properties as proposed 
long ago in Ref.s [1]. This energy loss can be quantified by the suppression of observed 
hadron spectra at high transverse momenta pt, namely Raa(pt) [2]. Altough the ob- 
servation of the jet suppression cannot be questionated there are several fundamental 
questions that still remain open. To investigate them we have constructed a model to 
study two observables beyond the Raa {.Pt ) • One is the elliptic flow v^ipr), that gives a 
measure of the angular dependence of quenching, and the other is Raa{q) / RAA(g) that 
determines the flavor dependence of the suppression. We suggest that the study of the 
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correlation between «2(j»t) and Raa{i) I Raa{9) carry information on the temperature 
dependence of the quenching and on the mechanism of parton flavor conversion. We find 
that an energy loss that increase as T — > T c , the q g in- medium conversion and an 
expansion-cooling of the fireball according to a lattice QCD EoS improve the agreement 
with the experimental data. Even if the T dependence that has to be considered in the 
present modeling including only path-lenth energy loss appear to be too extreme. 

2. Modelling the jet quenching 

In the model the density profile of the bulk is given by the standard Glauber model 
while the hard parton distributions in momenta space are calculated in the next-to- 
lcading- order (NLO) pQCD scheme. For the hadronization the Albino-Kramer-Kniehl 
(AKK) fragmentation functions have been employed. For further details sec Rcf. [3]. 
With regard to the jets energy loss we have employed various schemes, however to make 
a connection to the large amount of effort to evaluate gluon radiation in a pQCD frame, 
we have used also the Gyulassy-Levai-Vitev (GLV) formula at first order in the opacity 
expansion [4, 5]: 



where Cr is the Casimir factor equal to 4/3 for quarks and 3 for gluons, a s is the 
strong coupling, p(x, y, z) is the local density, r the proper time, E is the energy of the 
jet and fi = gT is the screening mass. There are corrections to Eq.(l) coming from higher 
order that can be approximately accounted for by a rescaling Z factor of the energy loss. 
However this is not really relevant for the objectives of the present work because we 
will renormalize the energy loss in order to have the measured amount of suppression 
Raa(pt) for central collisions. Usually in the GLV, as well as in other approaches, the 
temperature evolution of the strong coupling a s is discarded. We will consider the impact 
of such a dependence to understand the amount of T dependence coming simply from 
the asymptotic freedom. In the right panel of Fig.l we show by dot-dashed and dashed 
lines the temperature dependence of the energy loss for the GLV with a dependence of 
the coupling(GLV-a s (T)) and with a constant coupling a s = 0.27 (GLVc). Furthermore 
we show two other opposite cases for AE/ At: the thick line that shifts the energy loss 
to lower temperature (hence low density p or entropy density s) as suggested in [6, 7] 
and the other (thin line) that gives a dominance of quenching at high T, considered here 
just for comparison respect to the opposite case. We have applied our modelling of the 
jet quenching to Au+Au collisions at 200 AGeV and in the left panel of Fig.l we can 
see that, once the Raa{pt) at — 5% is fixed, the dependence on centrality is correctly 
predicted with a GLV formula for both constant and T-dependent a s . Hence looking at 
Raa one is not able to clearly discriminate the temperature dependence of the quenching 
and not even the details of the density profile [3] . 

3. Angular and Flavor dependence of the Quenching 

Generally, jet quenching modeling has not been able to simultaneously describe Raa 
and the elliptic flow v^. In particular, experimental data relative to v-i are considerably 



(1) 




1 

0.9 

_ 0.8 
> 

Q) fl7 

° 0.6 

v 

t±" 0.5 

v 

£ 0.4 
i0.3 



cc 



0.2 
0.1 






100 



200 

N „ 

part 



300 



400 



> 
O 

< 

LU 

< 



1.5 





GLV (a=00St) 

■ — GLV ce(T) J 

Energy loss at high T - 

Energy loss at low T I 


\ 

's 

: \ 

i . — r 


A *». \ - 



0.15 



0.2 0.25 
Temperature [GeV] 



0.3 



Fig. 1. - Left panel: Nuclear Modification Factor as a function of the number of participant 
in Au + Au at 200 AGeV [8] . Right panel: Temperature dependence of the energy loss for a 
parton with transverse momentum pr equal to 10 GeV. Dashed and dot-dashed lines represent 
the GLV energy loss with constant and T dependent a s coupling (see text). The thick line is 
the case in which the energy loss takes place only closer to the phase transition and the thin 
line represents an opposite case in which the energy loss take place only at high temperature T. 



larger than theoretical prediction. We have explored the relation between the temper- 
ature dependence of quenching and the value of the elliptic flow and in the left panel 
of Fig. 2 we can see that even if the amount of total quenching has been fixed to the 
experimental value of Raa{pt) the amount of elliptic flow is strongly dependent on the 
temperature dependence of the Ei oss . 

This correlation is due to the fact that at variance with Raa(pt) the i^Cpt) has 
a longer formation time because the jets have to explore the shape of the fireball to 
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Fig. 2. - Left panel: Elliptic Flow for pions (b=7.5 fm) coming from quark and gluon fragmen- 
tation for the different T dependence of the energy loss as shown in the right panel of Fig. 1. The 
shaded area shows the experimental data[8, 9]. Right panel: Ratio of quark to gluon Raa for the 
different T dependence of the energy loss as shown in the right panel of Fig.l. The shaded area 
approximatively shows the value expected for the ratio according to experimental observations 
and using the AKK fragmentation function 
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realize its asymmetry. Therefore it seems to be quite likely that experiments are telling 
us that quenching does not take place mainly at the very early time (temperature) of the 
collision but mainly at later times close to the phase transition. This would mean that 
the quenching is not proportional to the density (or entropy density) but a decreasing 
function of it with a maximum at T ~ T c . This is what is essentially discussed also in 
Ref.s[6, 7] where however it was implicitely assumed that the amount of quenching of 
quarks and gluons are equal among them and to the hadronic one. Here we have modified 
such assumptions showing that the temperature (or entropy density) dependence of Ei oss 
modifies not only the i>2(pt) but also the relative amount of quenching of quarks and 
gluons. 

3T. Quark to Gluon Modification Factor. - Due to its SU(3) Lie algebra the energy 
loss of gluons is 9/4 larger than the quark one. For this reason sometimes it is assumed 
that the ratio between the quark and gluon suppression (Raa(q) / Raa{9)) is equal to 
9/4. From this one would think that the (anti-)protons are more suppressed respect to 
pions because they come more from gluon fragmentation than from quarks fragmentation 
respect to pions, at least according to the AKK fragmentation function we employ. The 
data at RHIC however have shown that even outside the region where coalescence should 
be dominant [10, 11] the protons and the antiprotons appear to be less suppressed than 
the pions and p° [12, 13]. Again we can see that going beyond the simple amount of 
quenching given by Raa{pt) both the azimuthal dependence and the flavor dependence 
of the quenching appear to be in disagreement with the data. We call this open issues 
the "azimuthal" and the "flavor" puzzle respectively. We will show that even if Raa for 
central collisions is fixed to be ~ 0.2 the Raa{q) / Raa{9) is significantly affected by the 
temperature dependence of Ei oss . 

In the right panel of Fig. 2 we show the ratio of the Raa{<i) / Raa{9) for four different 
temperature dependences of the energy loss Ei oss , as in Fig.l (right). We can see that the 
standard GLVc energy loss does not give the expected ratio 9/4 for Raa{q) / 'Raa{9) but 
a lower value, around 1.8, which represents already a non negligible deviation from 2.25. 
We can however see that if the energy loss would be strongly T dependent and dominant 
in the T ~ T c region Raa(q) / Raa(9) can increase up to about 2.2 on the contrary, if 
it is dominant in the high temperature region (thin solid line) the Raa{q) / Raa{9) can 
become as small as 1.5. 

To understand this behavior is useful to consider the left panel of Fig. 3 where the trans- 
verse momentum distribution of initial parton and those obtained if all partons lose the 
same amount of energy (3 GeV and 4 GeV in the dotted and dot-dashed line, respec- 
tively) are shown. The effect of the quenching in this oversimplified case is to shift the 
spectra by a quantity equal to the amount of energy loss in a way indicated by the black 
arrow. Because of the rapid falling distribution the spectra after quenchig are in these 
cases one order of magnitude smaller respect to the initial one. Therefore the 10% of the 
spectrum without quenching, indicated by the thin line, is comparable to the spectrum 
in the case of Ei oss = 3-4 GeV . This means that if there are particles that lose a very 
small quantity of energy, like in the case of quenching at high temperature, they strongly 
influenced the final spectra for both quarks and gluons and damping the difference be- 
tween the Raa of quarks and gluons. For energy loss dominated by low temperature 
all particles lose energy and this increases the difference between the respective Raa ■ A 
similar effect could come not only from the T-dependence of Ei oss but from a core-corona 
effect [14]. For further explanations see Ref. [3]. 
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Fig. 3. - Left panel: Spectra for partons that lose a fixed amount of energy (see text). Right 
panel: Correlation between the Raa(q) / RAA(g) and the elliptic flow of pions with transverse 
momentum in the range 6 < pr{GeV) < fO. The squares refer to calculation without jet 
conversion while the circle are the values obtained including jet conversion with a K c = 6 
factor. 



4. — Correlation between Raa(q) / Raa(q) an d elliptic flow 

We have seen (Fig. 2, left and right panel) that an energy loss predominant at low 
T move toward observed data but move also the ratio (Raa (q)/Raa (<?) away from 
experimental indications. This is evident if one looks at the upper symbols of Fig. 3 
(right) where {Raa{q) / Raa{q) vs V2) is shown. To solve the "flavor puzzle" inelastic 
collisions that cause a change of the flavor has been invoked [15, 16, 17]. Such a process 
would at the end produce a net conversion of quarks into gluons. Hence a decrease of 
gluon suppression respect to the original suppression and an increase of the quark one. 
In Ref. [15] it has been calculated the conversion rate of a quark jet to a gluon jet and 
vice versa due to two-body scatterings. An enhancement factor K c = 4 — 6 that accounts 
for non-perturbative effect is needed to produce a nearly equal suppression of quarks 
and gluons. We have included such a mechanism in our model. The results are the lower 
symbols in the right panel of Fig 3. We can see that the q O g conversion does not affect 
the i>2 and allows to get closer to the experimental observed value, i.e. a V2 ~ 0.1 and 
an Raa{<i) / R.AA{g) < 1 (to account for the Raa(p + P) > Raa(k + + tt~) with AKK 
fragmentation function) . 

5. Impact of the Equation of State 

As a last point, we have observed that if the quenching is predominant near the 
phase transition the question of the corret equation of state arises. In fact the free gas 
approximation is no longer a reasonable approximation just close to T c . We have made 
some explorative studies on the impact that a more correct equation of state (EoS)can 
have on the correlation between Raa(q) / Raa{9) and v^ipr)- Making a fit to the lattice 
QCD data [18] we have obtained the following relation between density and temperature 



F. SCARDINA, ETC. 



0.4 
0.35 
0.3 
0.25 
u 0.2 
0.15 
0.1 
0.05 




GLVc 
GLV - a(T) 
E |oss at high T 
E. at low T 



WS-Glauber 



10 12 

p T [GeV] 



0.8 



0.6 



0.4 



0.2 



• b=9 fm 
■ b=7 fm 

♦ b=5 fm 
A b=3 fm 



4- 



0.5 



X 



♦ • 



4jrn/s 



QGP" 



1.5 
P T [GeV] 



2.5 



Fig. 4. - Left panel: Ratio between the elliptic flow and the eccentricity V2/e as a function of the 
transverse momentum pr of pions produced in collision with impact parameter equal to 7.5 fm 
for different temperature dependence of the energy loss, as shown in the right panel of Fig. 1. 
Right panel: «2/e as a transverse momentum function estimated in a parton cascade approach 
Ref. [20] at different impact parameters 



where /3(T) = 1/3 - a(T c /T) n with T > T c , a = 0.15 and n = 1.89 and of course 
for T >> T c one gets f$ <~ 1/3. In order to estimate the impact of this correction we 
have performed a simulation for the AEi oss (T) behavior represented by the thick solid 
line in the left panel of Fig.l which is similar to the delayed energy loss proposed by 
Pantucv [6] as a solution for the observed large elliptic flow. Wc consider only this case 
because it is of course the one that is much more affected by the modification implied 
by Eq.(2). The results are given by open symbols in the right panel of Fig. 3. Respect to 
the free gas expansion cooling the system spends more time at T ~ T c . This reflects in 
a further enhancement of both v-i (pt) and emcency of the q O g conversion moving the 
two observables closer to experimental data (shaded area in Fig. 3 (right)) 

As a last point we want to mention the issue of the efficiency of conversion of initial 
spatial asimmctry e into a i>2, namely the i>2/e. In the left panel of Fig. 4 we show 
the different values of v^j t that one obtains through the simple pathdength mechanism 
studied with our modeling, noticing that it is at maximum V2/e < 0.25. In the right 
panel of Fig. 4 we show the i^/e coming from elastic scattering in a cascade approach 
for a fluid at finite shear viscosity to entropy density Aitr]/ s = 1 [20]. Even if at slightly 
lower px the v-ij e in this case is about a factor of two larger. This seems to indicate that 
a proper treatment of elastic energy loss can give an important constribution at least at 
intermediate pr ~ 4 — 6 GeV where presently the experimental data on v-i are available. 



6. Conclusions 

We have pointed out the impact of peculiar temperature dependences of the energy 
loss on the elliptic flow and on the ratio between the quark and gluon suppression and 
their correlation. Moreover we have spot the relevance that the EoS may have in case 
of Ei oss dominant in the T ~ T c region. In any case our study, although already reveal- 
ing several interesting indications, is mainly explorative and it can be considered as a 
benchmark. A more quantitative analysis should be performed with more sophisticated 
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models that include the energy loss fluctuations, realistic gain and loss processes, elas- 
tic energy loss and a more accurate description of the bulk. Furthermore it should be 
explored if mass dependent formation time can affected the correlation between and 
RaaWi I Raa[.9) [21]. Obviouslyit is also important to study how the longer lifetime and 
higher temperatures which will be reached at LHC energies could affect the observed 
correlations. 
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